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Abstract. Using the brightest clusters in the ROSAT All-Sky Survey we compiled 
an X-ray flux-limited sample of galaxy clusters. The clusters have been reanalysed 
using ROSAT PSPC pointed observations if possible. The gravitational mass has 
been determined individually for each cluster in a homogeneous way assuming 
hydrostatic equilibrium. The mass function has been derived. We present the 
preliminary results and a comparison to previous determinations. 



1. Introduction 

Distribution functions of physical parameters of galaxy clusters can place important 
constraints on cosmological scenarios. Comparison of the mass function with 
analytical or numerical calculations can yield, e. g., the amplitude of the initial density 
fluctuations. Comparison of the X-ray luminosity, gas temperature or mass function 
in different redshift bins can give information about the cluster evolution. 

Several authors have published an X-ray luminosity function, e. g., jlL ||, Q, 
||. Also a cluster temperature function has been determined, e.g., [7j. A cluster 
gas mass function has been given by Burns et al. (1996), for an optically selected 
cluster sample. A gravitational mass function has previously been determined by 
Bahcall k Cen (1993), @, Biviano et al. (1993), |h}, and Girardi et al. Q. Bahcall 
& Cen used the galaxy richness and velocity dispersion to relate to cluster masses from 
the optical side and a temperature-mass relation to convert the temperature function 
of Henry & Arnaud (1991), ||, to a mass function from the X-ray side. Biviano 
et al. and Girardi et al. used velocity dispersions for an optically selected sample to 
determine the mass function. 

Since there is now high quality X-ray data available for the sample selection - using 
the ROSAT All-Sky Survey (RASS) - and detailed cluster analysis - using ROSAT 
and ASCA pointed observations - we have derived for the first time the galaxy cluster 
gravitational mass function using individually determined X-ray masses. 

This article will be published in the Proceedings of the 19 th Texas Symposium 
on Relativistic Astrophysics, held in Paris (1998) , and is also available at: 



http://www.xray.mpe.mpg.de/^reiprich/act/publi.htm 
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2. The Sample 

Completeness of a cluster sample is essential for the construction of the mass function. 
We compiled the clusters from RASS-based cluster surveys of high completeness 
(REFLEX, NORAS Q) and compared also with other published catalogs @, 
p5| , ]l6| ], |L7| , To avoid the high absorption and the crowded stellar field in which 
clusters are hardly recognized in the galactic plane, only clusters with a galactic 
latitude |6| > 20.0 have been included. For the same reasons the area around the 
Magellanic Clouds has been excluded. In addition the Virgo cluster region has been 
excluded here. The sky coverage is 26,720 deg 2 . 

We reanalysed the clusters using mainly ROSAT PSPC pointed observations and 
determined the X-ray flux /x(0.1-2.4 keV). 63 clusters have a flux greater than or equal 
to our adopted flux limit /x lim (0.1-2.4 keV)= 2.0 • 10 -11 erg/s/cm 2 . We call this cluster 
sample HiFluGCS (the Highest X-ray Flux Galaxy Cluster Sample). The distribution 
of HiFluGCS in galactic coordinates can be seen in Fig. |l|. 




Figure 1. Distribution of HiFluGCS in galactic coordinates. Open Triangles 
indicate clusters lying in the direction of the galactic plane, these are not included 
in HiFluGCS. 



3. Data Reduction and Analysis 

We used mainly high exposure ROSAT PSPC pointed observations to determine the 
surface brightness profiles of the clusters, excluding obvious point sources. If no 
pointed PSPC observations were available in the archive or if clusters were too large 
for the field of view of the PSPC we used RASS data. To calculate the gas density 
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profile the standard /3-modcl Q , [^9) (equ. |l|) has been used. 

/ r 2N-|/3 

PgasW -Pgas(O) f 1 + (1) 

Sx(R) =S X (0) ^ + 72-] ( 2 ) 

Fitting the corresponding surface brightness formula (equ. ||) to the observed surface 
brightness profiles gives the parameters needed to derive the gas density profile. 
To check if the often detected central excess emission (central surface brightness 
of a cluster exceeding the fit value) biases the mass determination we also fitted a 
double /3-model of the form Sx = Sx ± + Sx 2 and calculated the gas mass profile by 
/9 gas = ^/pg asi + Pg aS2 • Comparison of the single and double /3-model gas masses shows 
good agreement. 

We compiled the values for the gas temperature from the literature, giving 
preference to temperatures measured by the ASCA satellite p0| , plf , P^L fj|§ - 
For clusters for which we did not find a published temperature we used the X-ray 
luminosity-temperature relation given by Markevitch (1998), §. 

Assuming hydrostatic equilibrium the gravitational masses for the clusters can 
be determined. Plugging equ. |l| into the hydrostatic equation and assuming the 
intracluster gas to be isothermal yields the gravitational mass profile 

Having aquired the gravitational mass profiles for the clusters it is now important 
to determine the radius at which to determine the cluster mass. Simulations by 
Evrard et al. (1996), |24j, have shown that the assumption of hydrostatic equilibrium 
is generally valid within a radius where the mean gravitational mass density is greater 
than or equal to 500 times the critical density p c — 4.7 • 10~ 30 gcm~ 3 , as long as 
clusters undergoing strong merger events are excluded. This radius we call rsoo- We 
calculated the gravitational mass at r5oo and also r^m which is usually referred to as 
the virial radius. Using these definitions of the outer radius instead of a fixed length 
also allows the uniform treatment of clusters of different size. Using rsoo also saves 
us from an extrapolation much beyond the significantly measured cluster emission in 
general. 



4. Results 

In Fig. U we show the X-ray luminosity function of HiFluGCS compared to 
luminosity functions of other cluster samples. There is good agreement between these 
determinations, if anything than HiFluGCS shows a marginally higher density. 

In Fig. U we show the gravitational mass function for HiFluGCS for different 
definitions of the outer radius. Also shown are the mass functions obtained by 
Bahcall & Cen (1993), §], and Girardi et al. (1998), for an outer radius of 

3 /ig Mpc. Comparing our r5oo mass function with that of Bahcall & Cen we find 
increasing discrepancy towards lower mass clusters up to a factor of 7-8. For the 
r2oo mass function this discrepancy becomes less for the lower mass clusters but a 
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Figure 2. Luminosity function of HiFluGCS compared to luminosity functions of 
other cluster samples. 
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Figure 3. Cumulative gravitational mass function of HiFluGCS for different 
definitions of the outer radius compared to previous determinations which used a 
fixed radius of 3 Mpc. 



■5 



discrepancy arises towards the high mass end. In order to be able to directly compare 
the mass function for the clusters in HiFluGCS with the two others shown in Fig. || 
we determined the gravitational mass also at a fixed radius of 3 h^Q Mpc. Apart from 
the highest mass bin we almost exactly reproduce the mass function determined by 
Bahcall & Cen in this way. The value given by Girardi et al. lying a factor 3-4 higher. 

5. Discussion 

Two major points are of concern when deriving the mass function: 

1) The sample completeness and 2) the reliability of the mass estimates. 

1) We compiled the clusters from RASS-based cluster surveys. These surveys are 
complete at the 90 % level. The incompleteness of these surveys is likely to be highest 
at fluxes close to their adopted flux limit, which is much lower (factor ~ 5) than the 
flux limit adopted for HiFluGCS. Additionally we checked further published X-ray 
cluster catalogs. Besides the fact that still some clusters need to be checked in more 
detail we conclude that HiFluGCS is essentially complete. 

2) For the first time we have determined a mass function with cluster masses de- 
termined individually and in a homogeneous way for each cluster using high quality 
X-ray data. However, galaxy clusters are generally not spherically symmetric and 
hydrostatic equilibrium may not always be reached. Simulations by Schindler (1996), 
p5[ , and Evrard et al. (1996), [Q, however, have shown independently that the de- 
termined and true mass do not differ dramatically (< 20 %) if extreme merger clusters 
are excluded. Clusters also may not always be isothermal. For instance Markevitch 
et al. (1998), p0| , find a general trend that the temperature decreases with increasing 
radius, with the effect that the assumption of isothermality leeds to an overestimation 
of the gravitational mass in the outer parts (~ 6 core radii) of clusters by ~ 30 %. 
However, also average relative temperature profiles of cluster samples have been found 
which are consistent with being isothermal, e.g., Irwin et al. (1999), p6| . 

6. Conclusions 

By reanalysing the brightest clusters of RASS-based galaxy cluster surveys we 
constructed a complete X-ray flux- limited sample of galaxy clusters (HiFluGCS), the 
sky coverage being roughly 2/3 of the entire sky. We determined global physical 
parameters for the clusters using mainly high exposure ROSAT PSPC pointed 
observations. The luminosity function for HiFluGCS agrees well with previous 
determinations. We determined the mass function for the first time by individually 
determining the gravitational mass of each cluster in a homogeneous way assuming 
hydrostatic equilibrium and isothermality. Comparison with previous determinations 
shows a strong discrepancy especially towards lower mass clusters which is mainly due 
to the definition of the outer radius. For comparison we also determined the cluster 
masses at a fixed radius of 3 h^Q Mpc and apart from the highest mass bin we almost 
exactly reproduce the mass function determined by Bahcall & Cen. However, we 
suggest to use as outer boundary a radius which depends on the mean gravitational 
mass density, e. g. rsoo, J*200) m or der to treat clusters of different size in a comparable 
way. 

Another consequence of the definition of the outer radius becomes visible when one 
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integrates the mass function to determine the mass density bound in galaxy clusters 
Pbound- We find Abound relative to the critical density p c to be 1.0% for clusters of 
masses 2.5 • 10 13 h^Q M and higher using r 500 . The fraction increases slightly to 1.6 % 
for clusters of masses 4.1 • 10 13 h^ M Q and higher when we use r 2 oo- We find a larger 
increase if we formally calculate the fraction for a fixed radius of 3 h^ Mpc, which 
is 3.8% for clusters of masses 1.3 • 10 14 h$Q M and higher. Despite these different 
results depending on the outer boundary, however, it is clear that only a small portion 
of the total mass in the Universe is bound in galaxy clusters as the largest collapsed 
entities, implying that most of the mass must be somewhere else. 
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